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The XENON100 experiment, situated in the Laboratori Nazionali del Gran Sasso, aims at the 
direct detection of dark matter in the form of weakly interacting massive particles (WIMPs), based 
on their interactions with xenon nuclei in an ultra low background dual-phase time projection 
chamber. This paper describes the general methods developed for the analysis of the XENON100 
data, focusing on the 100.9 live days science run from which results on spin-independent elastic and 
inelastic WIMP-nucleon cross-sections have already been reported. 
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I. INTRODUCTION 

The XENON100 experiment aims to directly de- 
tect cold, non-baryonic dark matter which accounts for 
the majority of the matter in the Universe [1] accord- 
ing to plentiful astronomical and cosmological evidence. 
Among the most promising dark matter candidates are 
weakly interacting massive particles (WIMPs), arising 
naturally in several models beyond the Standard Model 
of particle physics [2] . If WIMPs are the dark matter 
particles, they could be directly detected via scatter- 
ing off nuclei [3]. In this paper, the complete analysis 
of the low-energy XENON100 data is presented focusing 
on the 100.9 live days science run. The XENON100 col- 
laboration previously reported exclusion limits on spin- 



independent elastic [4] and inelastic [5] WIMP-nucleon 
scattering using this data set. 

Following this introduction, Section II introduces the 
XENON100 experiment and details the formalism of the 
measurement process. Section III describes the data set, 
reconstruction procedures and selection cuts used in the 
analysis. This is followed by estimates of nuclear and 
electronic recoil backgrounds with control samples from 
Monte Carlo simulations, calibration data and science 
data outside of the signal region. Finally, a description of 
the observed event population and the comparison with 
the expected background is presented. 



II. THE XENONIOO EXPERIMENT 
A. Instrument description 
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The XENONIOO detector is filled with a total of 161 kg 
of ultra pure liquid xenon (LXe) divided in two concen- 
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trie cylindrical volumes. The inner target volume is a 
two-phase (liquid/gas) time projection chamber (TPC) 
of 30.5 cm height and 15.3 cm radius containing a xenon 
mass of 62 kg. It is optically separated from the sur- 
rounding LXe veto, closed on the bottom by a cathode 
mesh and on the top by a gate and an anode mesh, which 
provide the homogeneous electric fields required for the 
operation of the TPC. A technique similar to the one of 
a diving bell was chosen to keep the liquid in the TPC 
at a precise level between the gate and the anode meshes 
and to raise the liquid level in the veto (outside of the 
bell) above the TPC. This enables full LXe coverage of 
the TPC for an efficient suppression of external radioac- 
tive backgrounds. The XENON100 instrument and its 
subcomponents are explained in detail in [6] and only a 
short summary is given here. 

A particle interaction in the LXe target creates both 
excited and ionized atoms, which combine with the sur- 
rounding atoms to form excimers [7] . De-excitation of 
these excimers leads to a prompt xenon scintillation sig- 
nal (SI), which is recorded by photomultiplier tubes 
(PMTs) placed below the target in the LXe and above 
in the gas phase. Due to the presence of an electric drift 
field of 530 V/cm, a large fraction of the ionization elec- 
trons is drifted away from the interaction site in the TPC 
with a drift velocity ~ 1.73mm//is. Electrons which 
escape recombination and are not trapped by impurities 
are extracted from the liquid into the gas phase by a 
strong extraction field of ~ 12kV/cm, and a light sig- 
nal (S2) is generated by proportional scintillation in the 
gas [8]. The S2 signal is detected by the same PMTs 
but is delayed by the drift time td, which is the time it 
takes the electrons to drift from the interaction site to 
the liquid/gas interface. 3-dimensional event vertex re- 
construction is achieved using td and Vd to reconstruct 
the z position (z = Vdtd) and the hit pattern on the 
PMTs in the gas phase to reconstruct the (x, y) position. 
The ratio S2/S1 is different for electronic recoil events 
from interactions with the atomic electrons (from 7 and 
j3 backgrounds), and for interactions with the nucleus 
itself (nuclear recoils from WIMPs or neutrons), and is 
used to discriminate the signal against background. 

The PMTs are 1 inch x 1 inch Hamamatsu R8520, se- 
lected for high quantum efficiency (up to 32%) and very 
low intrinsic radioactivity. 80 tubes are immersed in 
the LXe below the TPC to ensure high light collection. 
98 PMTs are placed in the gas phase above the target, ar- 
ranged in concentric rings with the outmost ring extend- 
ing beyond the TPC edge for improved (x,y) position 
reconstruction. The LXe layer surrounding the target 
is instrumented with 64 additional PMTs, observing the 
volumes above, below and on the sides of the TPC, and 
operating as an active LXe veto. 

All materials used for the construction of the detector 
were selected for low intrinsic radioactivity [9] , in order 
to minimize the electronic recoil background and neu- 
tron background from (a, n) and spontaneous fission re- 
actions. To suppress the external background, the detec- 



tor is installed inside a passive shield. From the inside to 
the outside this consists of 5 cm of oxygen free high con- 
ductivity copper, 20 cm of polyethylene, 5 cm of lead with 
a low 210 Pb concentration, and 15 cm of standard lead. 
Three sides and the top of the shield have 20-cm-thick 
water or polyethylene shielding to further lower the neu- 
tron background, and the entire shield installation sits 
on a 20 cm slab of polyethylene. The detector is installed 
underground in the Laboratori Nazionali del Gran Sasso 
(LNGS) in Italy, at an average depth of 3 600 m water 
equivalent, which effectively reduces the muon flux by a 
factor of 10 6 compared to the surface flux [10]. 



B. Measurement process and event rates 

This section describes the formalism of light and charge 
detection in XENON100. Upper case latin letters in SI 
and S2 are used for the actual measured quantities while 
lower case letters in si and s2 denote expectation values. 
P is a discrete probability and p is a probability density 
function (pdf). 



1. Generation of light and charge 

For a given energy deposit of an interaction type u 
(u = nr for nuclear recoil and u = ee for electronic re- 
coil) E u in the presence of a drift field of strength £, 
the combined probability P(N 7 ,N e \E U ,£) for the gen- 
eration of N 7 photons and N e electrons shows, in general, 
an anti-correlation between the number of free electrons 
and photons. This is due to charge recombination pro- 
cesses which can lead to additional scintillation light [7] . 
At low energy deposits, such as for nuclear recoils in the 
dark matter region of interest, however, the measurement 
uncertainties dominate the width of the observed proba- 
bility distributions [11]. Hence, P can be approximated 
by independent Poisson processes: 

P(N 7 ,N e \E U ,£) w Poi(iV 7 I n 7 ) Poi(N e \ n e ). (1) 

The average energies needed for creation of one informa- 
tion carrier (photon or free electron) are expressed by 
effective "W- values", which depend on the interaction 
type, the drift field, and, at low energies, also on the 
deposited energy. The field dependence can be factor- 
ized with functions S u for the reduction in light yield 
due to field quenching and T u for the loss of charge due 
to recombination. The expectation values n~/(E u ,£) and 
n e (E u ,£) can then be written as 



n 7 (E u ,£) = 
n e (E u ,£) = 



E., 



E,, 



W^{E U ,£) W^(E u ,£ = 0) 

E u ^ E u 
W e (E u ,£) ~ W e (E u ,£ rci -»oo) 



Su(£), (2) 
T u (£), 



(3) 



3 



where 5^(0) = 1 and T u (£ rc f — > oo) = 1. 
In xenon dark matter detectors, the energy calibration of 
nuclear recoils is accomplished by comparing the signals 
from known 7-ray lines to dedicated measurements of 
the functions W 7 and S u , or W e and T u , which differ for 
electronic and nuclear recoils. In order to minimize the 
systematic uncertainties of the cross-calibration result- 
ing from modelling the detector responses, a reference 
source is frequently used. Historically, this has been the 
122keV co line from 57 Co decay, where the keV cc repre- 
sents the electronic-equivalent recoil energy. The in situ 
measured light and charge yields at 122 keV oc can be used 
as fixed points to establish the energy scale at lower 7- 
ray energies and for nuclear recoils, using the ratios of 
W 1 {E U ) and W e (E u ) relative to this reference, and ap- 
plying the functions S u (£) and T u (£), respectively. Cur- 
rently, for nuclear recoils (E u = E nr ), Wj(E nl ) [12] has 
been measured to lower energies than W e (E nI ). Hence, 
XENON100 uses the primary scintillation light to estab- 
lish the energy scale. 



where fi(f) — C y is the measured light yield for 

a reference 7-ray line (in PE/keV cc ) at the given electric 
field and position [see Eqs. (2) and (5)], given as 



C y (E cc — E rc f, £, r) 



S ee (£)n(r) 



WJE ee = E Iei ,£ = 0)' 



(7) 



C e H is the relative scintillation yield of nuclear recoils 
with respect to the reference 7-ray line at zero field, which 
equals 



C eS {E ni ,£ = 0) 



W 7 (E cc = E Iof ,£ = 0) 
W 7 (E ni .,£ = 0) 



(8) 



S nr and See are the reductions in light yield due to field 
quenching for nuclear and electronic recoils, respectively. 

The data analysis is usually performed with the 
"spatially-corrected" measured signal 



2. Measurement 0} primary scintillation light 

The expectation value for the primary scintillation 
light signal SI on PMT i with a gain gi 1 after integrating 
the current pulse is given in units of electron charge 

sl?(r) =n 1 (E U) £)j i (r)r) i g i = n 7 (E„, £)/Zj(f) g t , (4) 

where 7i(r) is the probability for a photon created at po- 
sition r within the TPC to reach the photocathode of 
PMT i and r)i is the product of quantum and collection 
efficiencies for that PMT. The combined function Hi(f) is 
the light detection efficiency and is measured with 7-ray 
calibrations (see Sec. IIIC). The raw data processor con- 
verts the measured signal into units of photoelectrons, 
using the estimates of gains which include an amplifica- 
tion factor of 10 [6]. Since errors in the determination of 
the combined PMT and amplification gains are typically 
< 2% [6] , these will be neglected in the following. Hence 
the expectation value for the primary scintillation signal 
on PMT i in units of photoelectrons (PE) is 



cSl = Sl(r) 



sU(r) w n 7 (E u ,£) m(r). 



(5) 



If M is the number of PMTs in the TPC, the energy 
deposit for nuclear recoils E nl determines the expected 
total primary scintillation signal si as 



=E nl C y (E cc = E Ic f, £, r) 

X C eS (Enr,£ = 0)^-, 



fi(r)' 



(9) 



(6) 



where n is the estimate for the sum over all PMTs and 
(//) is the spatial average. For such spatially-corrected 
signal, csl corresponds to the detector- averaged signal 
expectation value 



csl w n y (E ni ,£)(fi) 

= E ni (C y (E ieU £)) C cB (E ni ,£ = 0) §4S> ( 10 ) 

where (£ y ) is the detector-averaged light yield. 

Assuming a Poisson-distributed number of generated 
photons N 7 and a binomially distributed number of gen- 
erated photoelectrons for each PMT N pe i, the pdf is de- 
scribed by 

psi,i(Sli\n-y(E u ,£ )) dSl t = ^ y^P P mt,i(Slj | N peji ) dSlj 

N pCli N,, 

x Binom(7V peji | N 7 , ^(r)) 
x Poi(iV 7 \n 7 (E u ,£)) 

= Ppmt,i(Slj I Npe t i) 

x Poi(N P e ti I n 7 ,/ii(r))dSlt, 

(11) 

where p pm t,i(- ■ ■) is the response of PMT i, approxi- 
mated by a Gaussian with mean value N pe ^ and width 
cpMT, a/ N pe ^. The detected photoelectron width cpmT; 
is determined from PMT calibrations. 

Taking into account that SI = YliLi Sl< the pdf for 
the total light signal SI can be calculated from Eq. (11) 
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as 

Psi (Sl\n^E u ,£))dSl = 




(5(S1 -J2 Slj) dSli- ■ -rfSl M dSl. 

(12) 

This formula can be evaluated by computer simulations 
using individual response tables for each PMT. When 
considering only one spatially-averaged detector response 
where all PMTs are described by the same average re- 
sponse, Eq. (12) simplifies to 

PcS1 (cSl\n 7 (E u ,£))dcSl «^p pmt (cSl|iVp e ) 

N pe 

x Poi (iVp e I (/x)n 7 ) dcSl, 
(13) 

where N pe = J2i N pe ,i is the sum over all photo- 
electrons released on the PMT photocathodes. N pe 
is Poisson-distributed with the spatially-averaged mean 
value csl=n 7 (^i). p pm t(cSl|A^ e ) is the average PMT re- 
sponse function, approximated by a Gaussian with mean 
value N pe and width (JpMT\/N pe . The average singlc- 
photoelectron resolution of XENON100 PMTs is esti- 
mated to be <r pmt = 0.5 PE. 

3. Measurement of charge 



a narrow gas gap (h g ~ 2.5mm) , (3i(x,y) can be con- 
sidered as a function of (x,y) only. Gamma calibration 
lines can be used to measure the product 5i — kY (5^. 
Currently, only the sum over the PMTs is measured, re- 
sulting in an estimate S(x, y) = J2i $i( x ) Tj)- 

The analysis is usually performed with the corrected 
measured signal: 



cS2 = S2(r)e td/T = -f^-. (15) 
d[x,y) 

Using Eqs. (3) and (14), the expected total secondary 
scintillation signal for nuclear recoils can be written as 

s2(r) = s2i(r) = E nl Q y (E nl ) e -*-/ T « S(x, y), (16) 

i 

where Q y = T nr (£)/W e (E nr , £ rc f) is the measured charge 
yield of nuclear recoils (in e~/keV nr ) at the given electric 
field. The pdf ps2,i(S2i) can be described as 

Ps2,i(S2i|n e (i; M ,£ )) dS2i = ^ Ppmt,i(S2i | N pe:i ) 

x Poi(iV p(M | n e ,Si) dS2 i: 

(17) 

where p pmt ,i is the same response of PMT i as in Eq. (11). 
The pdf for the total proportional scintillation light sig- 
nal S2, denoted by ps2{S2\n e (E u , £)), can be evaluated 
analogously to Eq. (12). 



The ionization electrons produced at an interaction 
point drift through the liquid, where losses due to at- 
tachment to electronegative impurities with characteris- 
tic time r e occur (see Sec. IIIC). Electrons reaching the 
liquid surface are extracted into the gas phase with a 
yield k that depends on the extraction field £ gas . The 
same field is responsible for the proportional scintillation 
light signal S2 [8] resulting in an amplification factor Y 
from scattering of accelerated electrons with gas atoms. 
The expectation value for the secondary scintillation light 
signal on PMT i in units of PE is described by 



S 2 4 (r0 « n e {E u ,£) e'*'^' K (£ sas ) Y \^y,h^j fc(x, y) Vl . 

(14) 

Y is also called secondary scintillation gain and it de- 
pends on the ratio of £ gas to the gas density p, and on 
the size of the gas gap h g . Due to mesh warping or to an 
inclined liquid level, £ gaa and h g can be (x, y) position de- 
pendent. fii(x,y) is the probability for a photon created 
at the position (x,y) in the gas gap to reach the photo- 
cathode of the PMT i. Since the S2 signal is created in 



4- Event rate calculation 



The measured differential nuclear recoil rate 
d 2 R/dSldS2 for a given WIMP-nucleus scattering 



rate dR/dE nr is computed as 



' /_/ ' e(Sl,S2) f ^ p (Sl,S2\E nr )dE n 



dSldS2 v ' ' J dE, 
«ei, B (Sl)e 2> ,(S2) 
dR 



dE„ 



Psi {Sl\E m )p S2 (S2\E nr )dE nl ., (18) 



where ei jS (Sl) and £2, S (S2) are the SI and S2 signal de- 
tection efficiencies for the applied selection criteria on the 
data. These are a combination of the trigger threshold, 
event search algorithm and event selection cuts, where 
cuts can be applied on both spatially corrected and un- 
corrected SI and S2 signals. 

As already mentioned, so far the relation between E nT 
and SI signal has been measured more precisely and 
down to lower energies (n 7 (-E nr ,£) via £ e g(E nl )) than 
the relation between E nr and S2 signal (n e (E ny: ,£) via 
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Qy(Enr); see Fig. 5 of reference [13]). Therefore, the anal- 
ysis is done using only the measured differential rate ex- 
pressed in the cSl signal as 



dR 
dcSl 



'ei(cSl) 



dR 
dE nr 



p cS1 (cSl|£ m .) 



e 2 , s (S2)p S2 (S2|£; nr )dS2 dE n 



e2,E(E nr occsl) 

= £l(cSl) / -^e 2 . E (E m .)j5cSi(cSl \E nr )dE nr , 
J dE nr 

(19) 

where for most of the selection criteria the efficiency is es- 
timated in the spatially-corrected SI signal ei(cSl). The 
exception is the S2 threshold cut (S2 m i n ), which directly 
influences the corresponding SI signal efficiency via en- 
ergy sharing at the level of generated n 1 and n e before the 
Poisson fluctuations and the measurement uncertainties. 
Hence its acceptance e 2j E is calculated as a function of the 
detector-averaged expectation value csl (see Sec. Ill E 2), 
which is related to nuclear recoil energy E nr via Eq. (10). 
For notation simplicity, in the rest of the paper we use SI 
and S2 instead of the corrected quantities cSl and cS2, 
unless explicitly specified. 

III. DATA ANALYSIS 



anode had 9 trips due to the anode current exceeding a 
predefined threshold. Data from 20 minutes before the 
trip until 20 minutes after restoration of 4.5 kV were re- 
moved from the analysis, and the livetime was corrected 
accordingly. About 2% of the data was removed due to 
variations in the detector's operation parameters. 
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FIG. 1: Pressure and temperature of the XENON100 detec- 
tor during the science run. The dotted lines represent the 
maximum allowed variations and data periods falling outside 
are removed from the run selection. 



This section presents the complete analysis associated 
with the 100.9 live days dark matter search run acquired 
in 2010. The data set, analysis methods, event selec- 
tion and acceptances, and the background predictions 
are described in detail. This analysis has been used 
to obtain a result which has been already interpreted 
in terms of spin-independent elastic [4] and inelastic [5] 
WIMP-nucleon interactions. 



A. Data set 

The data used in this analysis were selected from pe- 
riods with stable detector operating conditions. Periods 
with xenon pressure and temperature values being more 
than 5 sigma away from the average value were excluded 
from the analysis. After this selection, other parame- 
ters like LXe level, cryostat vacuum pressure, purifica- 
tion flow and pulse tube refrigerator temperature were 
found to be stable during the whole science run. Fig- 
ure 1 shows the time evolution of the xenon pressure and 
temperature over the entire length of the run, including 
correlated fluctuations. A 5 sigma variation corresponds 
to ±0.04 K (0.04%) and ± 0.005 atm (0.24%) in temper- 
ature and pressure, respectively. 

The high voltages biasing the cathode at — 16 kV and 
anode at +4.5 kV were continuously monitored. The 
cathode was stable over the complete period while the 



The radon concentration in the XENON room and 
inside the shield cavity were measured using dedicated 
radon monitors 1 . The radon concentration in the 
XENON room was about 350Bq/m 3 (see Fig. 7 in [14]) 
and found to be stable for the whole period. The vol- 
ume inside the shield cavity was constantly flushed with 
boil-off nitrogen gas, resulting in a radon concentration 
< lBq/m 3 throughout the run, at the radon monitor 
detection limit and consistent with zero. The measured 
electronic recoil background rate was found to be stable 
over the whole time in various energy intervals. A pe- 
riod of 18 live days was removed from the data due to an 
increased level of electronic pick-up noise. 

The PMT gains were monitored with weekly LED cali- 
brations and found to be stable within 2% which is given 
by the precision of the measurement. Weekly 137 Cs cali- 
bration data (7-ray line at 662 keV ee ) were used to study 
the evolution of the light yield, charge yield, LXe purity, 
as well as the width of S2 pulses. The data was stable 
during the whole run and the LXe purity was continu- 
ously increasing (see Fig. 3). 

The science run started on January 13th, 2010 and 
ended June 8th, 2010. Taking into account time periods 
removed due to the above-mentioned instabilities, anode 



RAD7 from Durridge. 
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trips, high electronic noise levels and DAQ dead time 
(about 10% in this run), the total live time used for this 
analysis is 100.9 live days. 



B. Data acquisition and calibration 

For low-energy signals, the detector was triggered on 
S2 pulses using the analog sum of 68 PMTs in the in- 
ner part of the top array and 16 PMTs in the center of 
the bottom array. In order to trigger on very small S2 
pulses with typical 1 /is spread, the trigger signal was 
amplified, integrated and shaped with a spectroscopy 
amplifier, passed through a low threshold discriminator, 
and distributed simultaneously to the ADCs. The trig- 
ger threshold was measured directly using a dedicated 
background data. For this measurement, the discrim- 
inated trigger signal was split using a Fan-In/Fan-Out 
module, with one output distributed to the ADCs to 
trigger and the other one fed into an ADC channel to 
be recorded. The major S2 peak, which triggers the data 
acquisition, is usually followed by small S2 peaks (see 
Fig. 17). These mainly consist of single to few electrons 
events emitted in the gas phase due to photoionization 
of impurities in liquid xenon by the prompt or propor- 
tional scintillation photons [15]. Any additional trigger 
raised by the discriminator for any of the secondary S2 
peaks was recorded in the trigger channel. The ratio be- 
tween the S2 spectrum generating the signal in the trigger 
channel and the full S2 input spectrum gives the trigger 
function. Figure 2 shows the result of the trigger roll-off 
measurement. 
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FIG. 2: Trigger probability as a function of the S2 signal size. 

The background distribution due to electronic recoils 
is determined using the flat low-energy Compton contin- 
uum spectrum (see Sec. IIIF). A total of 5.8 live days of 
60 Co calibration data were acquired, spread in time over 
the run. The source was placed at three different posi- 
tions around the detector to evenly cover the sensitive 
LXe volume of the TPC. 

The nuclear recoil (signal) region is defined using elas- 
tic neutron interactions from 2.9 live days of 241 AmBe 
data, acquired right before the dark matter run. Most 
of the data quality cuts and their acceptance have been 
defined on these data (see Sec. HIE). 

The irradiation of LXe with neutrons also gives 7-ray 
lines at 40keV ce and 80keV ec from 129 Xe and 131 Xe, re- 



spectively. The xenon isotopes are also activated lead- 
ing to delayed de-excitation lines mainly at 164keV cc 
(t 1/2 = 11.8 d) and 236keV oc (t 1/2 = 8.9 d) from 131m Xe 
and 129m Xe, respectively. Some of these lines are used to 
define position corrections for the SI and S2 signals. 

The energy scale for nuclear recoils (E nr ) is inferred 
from the cSl signal via Eqs. (9) and (10). The scintil- 
lation efficiency C c s{E m ) of nuclear recoils relative to 
122keV co is parametrized using all existing direct mea- 
surements as shown in [4], The scintillation quench- 
ing factors due to the applied electric field for elec- 
tronic recoils See = 0.58 and nuclear recoils S nl = 0.95 
are taken from [16]. The light yield at 122keV co is 
(£ y ) = (2.20 ± 0.09) PE/keVcc and is interpolated using 
the light yields from the 40keV oc , 80keV cc , 164keV ee 
and 662keV oc lines [6]. 

The science data in the nuclear recoil region was 
blinded and therefore not accessible until the analysis was 
finalized. The lower 90% quantile of the S2/S1 discrim- 
ination parameter for 60 Co electronic recoils determined 
the blinding cut for energies Sl< 160 PE. The data out- 
side of the blinded region was also used for the estimation 
of the electronic recoil background in the WIMP signal 
region (see Sec. IIIF), for monitoring of the rate stability 
as well as for the estimation of the acceptances of certain 
data quality cuts. 



C. Data processing 

The digitized waveform from each of the 242 PMTs 
used in XENON100 are recorded. The trigger is located 
in the middle of the waveform, 400 jus long, and its global 
time is saved for each event. The raw data processor is 
based on the ROOT analysis toolkit [17]. It utilizes the 
difference in the SI and S2 peak width, where the former 
is of the order of a few tens of ns (see Fig. 6) mainly gov- 
erned by the LXe scintillation decay time constants [18] 
and the latter is of the order of a few /is determined by 
the electron cloud diffusion during the drift in the liq- 
uid and the drift time across the gas gap (see Fig. 9). 
For each peak candidate found by the processor, several 
peak properties are calculated and stored, for example 
pulse area (in PE), height (in V), width (in ADC sam- 
ples with 1 sample = 10 ns), etc. More details can be 
found in [6] [19]. 

For each identified S2, the (x,y) position is calculated 
using the hit pattern of the S2 photons in the top PMT 
array. Three independent algorithms were developed to 
reconstruct the position of each event. While an algo- 
rithm based on a neural network (NN) [20] is used for the 
analysis, the results from the two other algorithms (based 
on x 2 minimization and using support vector machines, 
SVM [21]) are used for quality cross checks, see [6]. 3 mm 
and 0.3 mm position resolutions (lcr) in (x,y) and in z 
are obtained, respectively. Because of the finite S2 signal 
width, only pulses which are more than 3 mm apart in z 
can be separated. 
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An energy deposit will produce a different number of 
detected SI photons depending on its position in the 
TPC. Primarily, this is due to the different paths that the 
light travels until a PMT is hit. A 2D axial-symmetric 
SI light collection map /Lt(r, z), determined from monoen- 
ergetic lines, is applied to each measured SI light pulse 
to obtain the position-corrected value [see Eq. (9)]. The 
maps derived using the 40 keV cc , 164 keV oe and 662 keV eo 
lines agree within 3%. 

Similarly, there are position-dependent effects in the 
detection of the S2 signal [see Eq. (15)]. A (x,y) light 
collection map S(x, y) is calculated, also from monoener- 
getic lines, using only the sum of S2 signals in the bottom 
PMT array (S2&) as it has a more homogeneous light col- 
lection efficiency. The maximum (x,y) correction in the 
48kg fiducial mass, selected for the analysis (see HIE 5), 
was 15%. An additional correction in z is necessary due 
to absorption of ionization electrons by electronegative 
impurities in LXe during their drift. The electron life- 
time r e , which is the time at which the total number of 
electrons produced is reduced by 1/e [see Eq. (14)], is 
used to quantify this effect [6] . The purity of the LXe 
was constantly increasing as monitored using the 137 Cs 
full absorption peak due to the continuous purification of 
the xenon gas by a hot getter. Figure 3 shows the time 
evolution of r e , which increased from 230 fis to 380 fis dur- 
ing this science run corresponding to a mean free path of 
40 cm to 66 cm. A linear fit is used to describe the data 
and the time-dependent S2 drift time correction was ap- 
plied accordingly. Since the maximum drift distance is 
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FIG. 3: Time evolution of the electron lifetime during the 
science run. A linear fit is used to describe the data. 

30 cm, the maximal drift time correction was smaller than 
a factor of 2. 



D. Analysis methods 

Two parallel analyses were performed to interpret 
these data in a spin-independent WIMP-nucleon inter- 
action framework as published in [4]. The energy region 
used for both analyses was (4-30) PE in SI (see HIE 2) 
corresponding to (8.4 — 44.6)keV nr - Before unblinding, 
it was decided to use the Profile Likelihood analysis 
method [22] as the primary interpretation method which 
does not employ a fixed discrimination in S2/S1 param- 
eter space. The calibration data sets were divided into 



12 bands along log 10 (S2b/Sl), in a way that for every SI 
bin the signal-like events are equally distributed between 
the bands (see Fig. 12). The number of bands was opti- 
mized based on the binning resolution of the discrimina- 
tion parameter log 10 (S2b/Sl) and the available amount 
of calibration data in each band. The total number of 
events from calibration data in the preselected energy 
region was about 4 300 for electronic recoils and about 
76 800 for nuclear recoils. While a large number of bands 
yields a result which depends less on the location of the 
event with respect to the band boundaries, the statistical 
uncertainty on the signal fraction in each band increases. 
It was verified that the results are robust to the number 
of bands, see also [22]. The background region consisting 
of Gaussian and anomalous contributions (see Sec. IIIF) 
is characterized using 60 Co calibration and non-blinded 
science data. The spectral shape of the neutron back- 
ground is taken from Monte Carlo simulations [23] . Each 
signal or background event has an associated probability 
to fall into a certain band. Based on the event distri- 
bution for signal and background on these bands, both 
the signal-plus-background and the background-only hy- 
potheses were tested regardless of the observed data. In 
the case of the exclusion of the signal hypothesis, the re- 
sulting limit on the WIMP-nucleon cross section is given 
at 90% confidence level (CL). 

A cut-based analysis was also performed in [4] to cross 
check the result. This analysis uses the optimum inter- 
val method [24] to derive a 90% CL limit on the WIMP- 
nucleon cross section. The observed number of events 
in the pre-defined log 10 (S2b/Sl) signal region (bench- 
mark WIMP search region) can be directly compared 
with the background prediction. For this analysis, af- 
ter all selection criteria applied, the benchmark region 
was constrained by the 99.75% electronic recoil rejection 
line defined using 60 Co data, the 3 a lower contour of the 
nuclear recoil distribution, the SI energy range and the 
threshold condition S2 > 300 PE (see Sec. HIE 5). 



E. Dark matter search event selection 

This section describes the criteria applied to the sci- 
ence data to select candidate events in the energy re- 
gion of interest, which fall into the following general cat- 
egories: data quality cuts, energy selection and threshold 
cut on S2, selection of single scatter events, consistency 
cuts, selection of the fiducial volume, and the signal re- 
gion selection for the cut-based analysis. 

The acceptance loss for each condition is determined 
using the fraction of events removed by this selection 
alone. For the case of independent selection cuts, this 
method gives a better estimate of the acceptance than 
using the fraction of the total events which pass the cut 
under study. The reason is that events which fail mul- 
tiple cuts are far less likely to be valid events. For the 
cases in which the selection cuts were found not to be 
independent, their combined acceptance is derived. Vi- 



sual inspection of a subset of both accepted and rejected 
event waveforms was used as a cross-check to validate 
the acceptance determination. The acceptance was eval- 
uated using the energy range and the 48 kg fiducial mass, 
selected for the analysis (see HIE 2 and HIE 5). Nuclear 
recoil data is used to determine the acceptance for most 
of the cuts, with the few exceptions where conditions of 
neutron calibration datasets were found not to be repre- 
sentative for the full science data taking conditions. 



1. Data quality 

Data quality cuts remove non-physical or noisy wave- 
forms. A two-fold PMT time-coincidence within a 20 ns 
time window is required for a valid SI, where each PMT 
hit must have a signal size larger than 0.35 PE. This effi- 
ciently rejects PMT dark current signals as they are un- 
likely to happen simultaneously in multiple PMTs. Ad- 
ditionally, the SI coincidence requirement was tightened 
for the pulses where the SI peak contains noisy PMT 
channels. If known noisy PMTs contribute to the SI co- 
incidence, then the coincidence requirement is increased 
by the number of noisy PMTs [4] to avoid identification 
of an SI candidate with only one good PMT. 

For the results presented in [4, 5] , the acceptance was 
computed via a Monte Carlo simulation, which takes into 
account light collection in the TPC, PMT quantum effi- 
ciencies, single photoelectron resolution, and the time de- 
pendence of the scintillation process. Using this method, 
the acceptance was estimated to be > 99% above 6 PE 
and dropping down to 97% at 4PE. However, the ac- 
ceptance of this cut can be evaluated using data instead 
of a Monte Carlo simulation. Events from nuclear re- 
coil calibration data with a coincident signal in the veto 
are used to compute the acceptance of the coincidence 
requirement, as they correspond to true physical interac- 
tions in the TPC. The acceptance is > 99% above 10 PE 
and drops down to 80% at 4PE (see Fig. 4). This new 
method gives a slightly lower acceptance compared to the 
previous one, but this has a negligible impact on the in- 
terpretation of the results in [4, 5]. Since the data-driven 
method gives a better estimate of the cut acceptance it 
is used in the upcoming data analysis. 

A population of events with an SI peak due to noise 
was found after unblinding in the region around and be- 
low the analysis threshold of 4PE. This population is 
due to electronic pick-up noise and not related to phys- 
ical events in the TPC. Since the original cuts were not 
sufficiently restrictive against electronic noise, two addi- 
tional cuts had to be introduced post-unblinding in order 
to remove the full noise population, including three noise 
events with SI signal leaking above 4PE. The SI coin- 
cidence requirement was modified by adding few more 
noisy and non-functional PMT channels. The latter are 
turned off but can still pick-up electronic noise. An ex- 
ample of a rejected noise event is shown in Fig. 5. The 
modification of the SI coincidence requirement had no 
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FIG. 4: Acceptance of some of the cuts as a function of SI: 
SI coincidence requirement using the published Monte-Carlo- 
derived (solid blue) and new data-driven (dotted blue) meth- 
ods, waveform quality cut using SI and S2 candidate peak 
area (red dash-dot-dotted), S2 single scattering (green dash- 
dotted), and S2 pulse width cut (magenta dashed). The ver- 
tical lines represent the energy window for the analysis, see 
Sect. HIE 2. 



relevant impact on the nuclear recoil acceptance (< 0.4% 
acceptance loss) as determined on 241 AmBe calibration 
data. 
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FIG. 5: PMT traces for a noise event found after unblinding. 
The SI candidate region is indicated by blue dashed lines. All 
PMT signals contributing to the SI candidate are periodic 
electronic pick-up noise peaks. The bottom panel shows a 
non functional PMT. 

In addition, a cut on the SI width was introduced. 
Figure 6 shows the histogram of the SI peak width, de- 
fined between the right and left peak boundaries at 10% 
of the height. The distribution has a regular structure 
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FIG. 6: Histogram of the SI width at 10% level for neutron 
data. The line represents the cut on this parameter, events 
left of the line are rejected to eliminate electronic noise. 
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due to the finite sampling frequency of 100 MHz. Very 
small widths correspond to electronic noise pulses with 
samples fluctuating around the baseline, shortening the 
extent of the identified peak due to negative excursions 
(see Fig. 5). Events left of the vertical line are rejected. 
The line was defined after visual inspection of waveforms 
from 241 AmBe and low energy science data. Waveforms 
of valid events passing all cuts are shown in Fig. 17. The 
SI width cut has an acceptance of 99.9% at 4PE in- 
creasing to >99.9% above 6 PE, as verified with nuclear 
recoils. 

The peak finding algorithm is tuned to correctly iden- 
tify low-energy interactions with the highest possible ac- 
ceptance. In waveforms of events with large energy de- 
positions in the TPC (e.g. through-going muons) or with 
HV-induced micro-discharges, features in the pulse shape 
can sometimes be mistaken by the algorithm as individ- 
ual small energy signals. The event selection requires 
that the integral of the largest SI and S2 signals is big- 
ger than the size of the integral of the remaining wave- 
form. In order to obtain a sample with DAQ conditions 
representative of the full science run, the acceptance of 
this cut was calculated using science data outside of the 
blinded region. The cut acceptance in the energy region 
of interest is about 98% (see Fig. 4). 

There are events where no valid SI peak is present in 
the waveform. Such events appear cither if the real SI 
is too small to be detected or is missed by the peak- 
finder, or if it is too close to the S2 peak to be resolved 
as an individual peak. The latter mostly happens for 
interactions at the very top layer of LXe or in the gas 
region. In these events, a random electronic noise or 
a coincident PMT dark current peak can be picked up 
as an SI candidate. Their very asymmetric S2 signal 
distribution on the top and the bottom PMT arrays can 
be used to identify and reject them. This selection was 
defined using the non-blinded part of the science data. 
The cut acceptance was determined to be ~ 99.6% in the 
energy range of (4 — 30) PE using neutron calibration 
data. 

Occasionally, S2-like events with signals very localized 
in the (x,y) plane and/or unusual light patterns occur. 
Two classes of such spurious events were found in some 
of the science data outside of the blinded region. In the 
first class, the events have the bulk of S2 light detected 
by a single PMT. As the maximal fraction of the S2 light 
seen by a single PMT has a mean value of about 20%, 
events are rejected when this fraction is more than 65%. 
In the second class, events have an unusually small frac- 
tion of S2 light detected by the top PMT array. Since 
the S2 signal is generated in the gas phase, the top PMT 
array typically sees about 1.3 times more proportional 
scintillation light than the bottom array. Hence events 
are removed when the amount of S2 light detected by 
the top array is less than 74% of the light in the bottom 
array. Both cuts were optimized using the non-blinded 
science data. Their acceptance was calculated using nu- 
clear recoil data and it is above 99.6%. 



2. Energy selection and S2 signal threshold 

As already mentioned, the SI energy region used for 
the analysis ranges from 4 to 30 PE. The selection ac- 
ceptances above 4 PE were found to be high and nearly 
energy independent, while no notable WIMP sensitivity 
was gained by increasing the energy window above 30 PE. 
This region corresponds to (8.4—44.6) keV nr (see Eq. 10). 
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FIG. 7: Acceptance of the S2 threshold cut as a function 
of csl, the SI signal before Poisson fluctuations. The sharp 
drop at 1 PE was introduced artificially in order to be more 
conservative at the lowest signals. This acceptance is used in 
the WIMP rate calculation before the SI signal is convoluted 
with a Poisson function to take statistical fluctuations into 
account. 

In order to stay well above the trigger threshold which 
starts to roll off around 280 PE (see Fig. 2) a valid event 
is required to have S2>300PE (uncorrected S2). The 
acceptance of such a condition affects the measured dif- 
ferential rate [see Eq. (19)] and is determined in the fol- 
lowing way. For a given interval of position-corrected 
cSl from 241 AmBe data, the position-corrected cS2 spec- 
trum is described in a functional form and extrapolated 
conservatively to the range below threshold. In order to 
minimize the extrapolation, only the top 7.5 cm of the 
detector are used for this analysis reducing the impact 
of charge losses due to the finite r e . From this data the 
acceptance £2(cSl) of the S2 > 300 PE condition is deter- 
mined numerically, taking into account the 48 kg fiducial 
mass, the increasing r e during science data taking, and 
the spatially uniform event distribution as expected from 
WIMP interactions. 

This acceptance, however, cannot be simply applied 
to the observed science data as the measurement-related 
fluctuations of the SI and S2 signals are independent 
from each other (see Sect. II B). As an example, a recoil 
from a 7 GeV/c 2 WIMP can enter the detection region of 
(8.4 — 44.6) keVnr only because of upward fluctuations of 
the SI scintillation signal. The average S2 signal of these 
events, however, is smaller than from WIMPs of higher 
masses, leading to a lower acceptance at the same SI sig- 
nal. The acceptance of the S2 threshold cut is therefore 
calculated as a function of the csl signal without Poisson 
fluctuations, which is proportional to the nuclear recoil 
energy deposition E nr [see Eqs. (6)-(10)]. Since the nu- 
clear recoil spectrum varies with the WIMP mass, this 
leads to a WIMP-mass-dependent acceptance (see Fig. 2 
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of [4]). 

In general, csl and therefore the acceptance E2,e(cs1) 
are experimentally not accessible. However, the latter 
can be determined indirectly from the previously derived 
e 2 (cSl) from 241 AmBe data. Using C e g(E nT ) from [4] and 
Eq. (10), a 241 AmBc Monte Carlo simulated energy spec- 
trum is converted into a csl signal spectrum. An initially 
estimated acceptance function e* 2 E (csl) is now applied to 
this spectrum and varied iteratively until the acceptance 
e^cSl) from the Monte Carlo dataset convoluted with a 
Poissonian distribution to account for the statistical fluc- 
tuations is equal to e 2 (cSl). Figure 7 shows the resulting 
acceptance £2,e(cs1). The acceptance is conservatively 
set to be below l PE. 



3. Selection of single scatter events 

WIMPs are expected to scatter only once in the TPC 
producing one SI and one S2 pulse in the waveform and 
no coincident signal in the veto volume. S2 pulses of a 
few PE size after the main S2 signal or small ionization 
signals from single electrons can bias the single scatter 
event selection as they can be mistaken as an S2 signal 
from a second scatter. It was observed that the size of 
these delayed signals is correlated with the size of the 
main S2 pulse. Thus, events with one S2 peak are se- 
lected by requiring all other S2 peaks to be smaller than 
a certain threshold value which was defined using nu- 
clear recoil calibration data. However, the acceptance 
was calculated using electronic recoil science data, since 
the low-energy region was dominated by beta particle in- 
teractions which represent an almost pure single scatter 
event sample (see Sec. IIIF2). The cut acceptance in 
the energy region of interest is > 95% (see Fig. 4) and is 
independent of the S2 signal size. 

Even in the case of an event with multiple scatters, 
only one SI peak is expected since the 10 ns sampling 
time does not allow to distinguish between them. Multi- 
ple SI signals in a waveform would come either from pile 
up (which is very unlikely given the low trigger rate of 
~ 1 Hz during science data taking), PMT dark current, 
electronic noise, or misidentified uncorrelated S2 ioniza- 
tion signals from single electrons. Events where a second 
SI peak is seen by at least two PMTs are therefore re- 
moved if it appears within the maximum possible drift 
time. However, if the S2/S1 ratio is too high for the sec- 
ond SI candidate, the event is not rejected. Since this 
cut is affected by the noise conditions during the mea- 
surement, the acceptance was determined from the non- 
blinded science data to be ~ 98.5% in the (4 - 30) PE 
energy range. 

Given the negligible rate of accidental coincidences, a 
signal in the active LXe veto which is in coincidence with 
the SI peak in the TPC must therefore be from a double 
scatter interaction. All events where a coincident signal 
in the veto is present with a size > 0.35 PE are rejected. 
The non-blinded science data is used to calculate the ac- 



ceptance of the cut, which is > 99.6%. 

4- Consistency cuts 

Events were selected based on the comparison of the 
reconstructed position estimators from different algo- 
rithms, the fit quality of the reconstructed position, and 
the comparison of the measured SI PMT pattern with 
the one expected from the reconstructed position. Inter- 
actions in the gaseous xenon or events with SI and S2 
signals not from the same physical interaction were also 
removed. 




FIG. 8: Top: Neural network Xr value versus the combined 
difference in the reconstructed position by the NN, SVM and 
X 2 — minimization reconstruction algorithm from the 241 AmBe 
calibration data. Bottom: example of a double scatter event 
from neutron data rejected by the position reconstruction 
cuts. The PMT pattern of the top array clearly shows the 
presence of two S2 signals centered roughly on PMT 32 and 
45, respectively. The black dot close to the center represents 
the reconstructed position. The color code of the legend repre- 
sents the measured S2 signal size (in PE) seen by each PMTs. 

Large differences in the event positions inferred by the 
three different algorithms (see Sec. Ill C) in general corre- 
spond to unresolved multiple scatter events as their inter- 
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actions are so close in z that the corresponding S2 signals 
cannot be separated. Events with the combined differ- 
ence, calculated as the square-root of the sum of squared 
differences between (x, y) positions reconstructed by NN, 
SVM and ^-minimization algorithms, larger than 7 mm 
are removed. A % 2 -value of the difference between the 
observed S2 PMT hit pattern and the one expected from 
the Monte Carlo simulation quantifies the quality of the 
reconstruction position. Hence, a cut on the reduced Xr 
value (x 2 divided by the number of PMTs in the upper 
array) of the NN algorithm is used to reject badly recon- 
structed events. Figure 8 (top) shows the Xr distribution 
of events calculated from an NN algorithm versus the 
combined difference in reconstructed position. Events 
above and on the right of the solid lines, are rejected. 
Figure 8 (bottom) shows the S2 PMT pattern of a dou- 
ble scatter event which is rejected by these conditions. 
The position is reconstructed wrongly (black dot) close 
to PMT 97, in between the two true positions. The ac- 
ceptance of these conditions is evaluated using 60 Co and 
241 AmBe calibration data. It is > 99.6% in the energy 
range of (4 - 30) PE. 

The width of an S2 pulse increases with the td due to 
the diffusion of the electron cloud during its propagation 
towards the gate grid. Typical width values, defined at 
10% of its peak height, range from ~ 1 /zs to ~ 2/zs for 
minimal and maximal (td = 176 fxs) drift times, respec- 
tively. The independent S2 width is compared to the td 
value and events are rejected when these quantities are 
inconsistent. Neutron calibration data are used for the 
cut definition. Low energy S2 events show larger spread 
due to low statistics of drifted ionization electrons, hence 
the cut is defined in an energy-dependent way. The cut 
acceptance is set to 90% (see Fig. 4). The top panel of 
Fig. 9 shows the distribution of the S2 width versus td for 
nuclear recoil events from neutrons for all energies (blue). 
The red dots are the events removed by the S2 width se- 
lection. The waveforms in the middle and bottom panels 
correspond to accepted typical candidates with td = 8 /is 
and td = 127 /is between the SI and S2 peak. Their un- 
corrected S2 sizes are ~1 000 PE and ~1 500 PE, respec- 
tively. 

Multiple scatters from electromagnetic interactions in 
the science data or 60 Co calibration data, where one 
interaction happens inside the TPC and at least one 
in a charge insensitive region, can produce anomalous 
electronic-recoil leakage events. Such regions include the 
LXe volume below the cathode or close to the field shap- 
ing rings near the inner TPC wall. These events, referred 
to as "gamma-X", are classified as valid single scatters 
because a single S2 is measured but the SI is the sum 
of both prompt SI signals. Having an S2/S1 ratio lower 
than that of true single scatters, the events can poten- 
tially leak into the WIMP search region. 

The PMT hit pattern of the SI signal is used to dis- 
criminate between these events and true single interac- 
tions. A likelihood parameter is computed as —2 log(Ap), 
where Xp is the ratio of the Poisson likelihood of the 
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FIG. 9: Top: Distribution of the S2 width (defined at 10 % 
of the peak height) versus drift time for nuclear recoils with 
energies up to 30 PE in SI (blue histogram). The red dots 
are the events removed by the S2 width selection. Middle: 
S2 waveform example with td =8 fis. Bottom: S2 waveform 
example with td = 127 /is. 



measured SI PMT pattern and the "standard" SI PMT 
pattern of single scatter events which happen at the same 
reconstructed (x, y, z) position. The standard SI PMT 
pattern map is acquired from the full absorption peak of 
137 Cs calibration data. Figure 10 shows the distribution 
of the likelihood parameter for neutron data as a func- 
tion of the uncorrected SI signal area. Events with a poor 
likelihood ratio (above the line) are rejected because they 
have an SI pattern that is inconsistent with the expecta- 
tion based on their position. The acceptance of this cut 
is determined using 60 Co and 241 AmBe calibration data 
and is ~ 97% with no energy dependence. 

Figure 11 shows an example of an event removed by 
the pattern likelihood selection which could be explained 
by a double scatter event with one interaction below the 
cathode. The top panel contains the pattern of the S2 
signal in the top PMT array. The position is well re- 
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constructed close to PMTs 72 and 87 (black dot). The 
bottom panel shows the pattern of the SI signal in the 
bottom PMT array. The reconstructed z position is very 
close to the bottom PMT array, about 2.5 cm above the 
cathode grid and therefore the major part of the SI light 
signal is seen around the (x,y) reconstructed position. 
In addition, there is a large portion of the SI signal seen 
right above PMT 113 with no charge detected at simi- 
lar (x, y) pointing to a possible second scatter below the 
cathode. 
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FIG. 10: Distribution of the likelihood parameter of the PMT 
pattern as a function of the uncorrected SI signal for nuclear 
recoils. The likelihood parameter is computed as — 21og(Ap), 
where Ap is the ratio of the Poisson likelihood of the measured 
SI PMT pattern and the "standard" SI PMT pattern at the 
same position (see text). The red line represents the cut line 
above which events are discarded. 



5. Signal/background discrimination and fiducial volume 

Figure 12 shows the electronic and nuclear recoil bands 
from 60 Co (blue) and 241 AmBe (cyan) data, respectively, 
using the discrimination parameter log 10 (S2b/Sl), flat- 
tened by subtracting the electronic recoil mean. This 
removes the energy dependence of the electronic recoil 
band allowing for an easier combination of data from dif- 
ferent energies. The border lines of the 12 bands used 
by the Profile Likelihood analysis to model signal and 
background distributions are shown in red in the chosen 
energy window. 

Compared to the Profile Likelihood analysis, the 
benchmark WIMP region for the cut-based analysis is 
additionally constrained by the S2/S1 electronic recoil 
rejection cut and the lower bound of the nuclear recoil 
band. The electronic recoil rejection level was defined 
using 60 Co calibration data which has an S2/S1 distri- 
bution is described well by a Gaussian function with an 
extra tail appearing about 3 a away from the mean value. 
It was set to 99.75% (see Fig. 12). The acceptance of this 
cut, calculated using nuclear recoils, is energy dependent 




FIG. 11: Example of an event removed by the pattern like- 
lihood selection. Top: The top array PMT pattern of the 
S2 signal with the reconstructed (x, y) position (black dot) . 
Bottom: The bottom array PMT pattern of the SI signal. 
This could be due to a double scatter event. One interaction 
happens below the cathode and close to the bottom PMT ar- 
ray resulting in a very localized SI above PMT 113 with no 
corresponding S2 signal at similar (x,y). The second inter- 
action happens in the sensitive volume but still close to the 
cathode grid. Its (x, y) position is well reconstructed using 
the S2 signal in the top array, and its SI signal is localized in 
several PMTs in the bottom array at similar (x,y) position, 
due to the vicinity of the interaction point. 



and ranges from ~ 29% to ~ 46% (see Fig. 2 of [4]). The 
241 AmBe data is used to constrain the signal region re- 
moving events which are more than 3 a away from the 
mean of the nuclear recoil band (see Fig. 12). 

An ellipsoid containing 48 kg of liquid xenon is used 
as the fiducial target mass. Figure 13 shows the shape 
of the fiducial volume cut (red) together with the ob- 
served event distribution from the science data without 
applying any electronic recoil discrimination cut. The 
yellow dashed-contour represent the borders of the TPC. 
Events falling outside of it are due to the uncertainty in 
the position reconstruction. The background during this 
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tector are excluded as the position corrections are the 
largest in those regions. A few milimeters near the grids 
are removed to avoid the increased background in this 
region. 

The log 10 (S2b/Sl) distribution shows a population of 
events at low energies (few PE) and high values of the dis- 
crimination parameter, log 10 (S2b/Sl)— ER mean > 0.5. 
These correspond to random small SI peaks that do not 
belong to the S2 signal. A dedicated cut removes these 
events which are far above the nuclear recoil band with 
an acceptance above 99.95%. 

6. Overall cut acceptance 



FIG. 12: Electronic (blue) and nuclear (cyan) recoil bands in 
flattened discrimination space. The black dashed lines, rep- 
resenting the analysis energy window (vertical) and the S2 
threshold cut (bottom left), are used in both analysis meth- 
ods. The 12 bands used in the Profile Likelihood analysis to 
model signal and background are shown in red. The bench- 
mark WIMP region for the cut-based analysis is further con- 
strained by the black dotted lines which are 99.75% electronic 
recoil rejection line (middle) and the lower 3 a contour of the 
neutron distribution (bottom right). 




FIG. 13: Distribution of all events in the TPC observed in the 
8.4 — 44.6keV nr energy range during 100.9 live days. With 
the exception of the electronic recoil discrimination cut, all 
cuts are used including the ones introduced post-unblinding 
to remove a population due to electronic noise. The fiducial 
volume containing 48 kg of liquid xenon (red) and the TPC 
borders (yellow) are also indicated. Events outside of the bor- 
ders are due to the uncertainty in the position reconstruction. 



Figure 14 shows the overall acceptance of all cuts used 
for the Profile Likelihood analysis in the (4 - 30) PE SI 
energy range and 48 kg fiducial mass. The solid line in- 
cludes the published SI coincidence requirement derived 
by a Monte Carlo simulation and the dotted line the im- 
proved acceptance using a data-driven method, see Fig. 4. 
The acceptance of the S2 threshold is shown in Fig. 7. 
The later has to be applied first to the WIMP spectrum 
without statistical fluctuations. Then the combination of 
all other cuts as shown in Fig. 14 is applied to the SI spec- 
trum after taking into account Poisson fluctuations. As 
already mentioned, for the cut-based analysis, the S2/S1 
electronic recoil rejection cut and the 3 a lower nuclear 
contour are additionally applied (see Fig. 2 of [4]). 
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FIG. 14: Measured acceptance of all quality cuts as a func- 
tion of the corrected SI signal using the published SI coinci- 
dence requirement derived by Monte Carlo (solid blue) and 
by a new data-driven method (dotted blue), see Fig. 4. The 
S2 threshold cut is treated separately (see Fig. 7). The two 
vertical dashed lines indicate the pre-defined energy window 
(4 - 30) PE =(8.4 - 44.6) keV nr . 



F. Background expectation 



science run was dominated by 85 Kr which is uniformly 
distributed in the TPC as is visible in the central region 
of Fig. 13 (see also Sec. IIIF2). For this reason, a smaller 
fiducial volume would not help to reduce the background 
further, and the largest possible volume was used instead. 
This choice was based mainly on the understanding of the 
spatial corrections in the detector. The edges of the de- 



This section describes the various background con- 
tributions and quantifies the number of expected back- 
ground events for the 100.9 live days exposure, the 48 kg 
fiducial mass and the benchmark region of the cut-based 
analysis defined in Sec. HID. The Profile Likelihood anal- 
ysis uses exactly the same data and assumptions for the 
background estimation. 
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1. Nuclear recoil background from neutrons 

The nuclear recoil background from muon-induced 
neutrons and from neutrons created in spontaneous fis- 
sion and (a,n) reactions in the detector materials are de- 
termined using a GEANT4-based [25] Monte Carlo sim- 
ulation [23] . Single scatter nuclear recoils in the LXe 
with no associated electronic recoil from inelastic reac- 
tions were selected assuming 3 mm event resolution. 

Cosmogenic muons were simulated using the MU- 
SIC [26] and MUSUN [27] packages. The muons and their 
daughter particles from the electromagnetic and hadronic 
showers were propagated through the rock and through 
the XENON100 shield and detector materials. A total of 
(0.08^0 04) muon-induced events were found. 

To obtain the neutron rate from natural radioactiv- 
ity, the a-activity from the uranium and thorium chains 
is considered. These are derived from the detector 
material's 7-activities measured with germanium spec- 
troscopy [9] . The neutron production rate in the de- 
tector and shield materials is computed using the a- 
activities from a modified SOURCES4A code [28, 29]. It 
calculates neutron production rates and energy spectra 
from (a,n) reactions and spontaneous fission. A total of 
(0.032 ± 0.006) neutron background events due to radio- 
genic background were found. 

Summing up all contributions, the total nuclear recoil 
background is (0.1lt°;° 8 ) events, dominated by cosmo- 
genic neutrons. The rather large errors are due to sys- 
tematic uncertainties in the neutron production rates in 
the different materials [23] . 



2. Electronic recoil background from 7 and j3 events 

During the 100.9 live days run, the background at low 
energies was dominated by intrinsic 85 Kr (beta with 
687keV co endpoint), introduced accidentally through a 
tiny air leak in the gas system before the start of the 
dark matter search. The 85 Kr concentration was deter- 
mined by comparing the measured energy spectrum to 
a detailed Monte Carlo simulation. The simulation in- 
cludes the intrinsic radioactivity of all detector compo- 
nents, based on the measured contamination of the vari- 
ous materials. This method is described in detail in [14]. 
Assuming a 85 Kr isotopic abundance of 10 -11 , a natural 
Kr concentration of (700 ± 100) ppt was obtained for the 
data presented here. In parallel, an analysis using de- 
layed coincidences was performed to determine the 85 Kr 
concentration with an independent method [30] . It uses 
the decay of 85 Kr to 85m Rb (beta with 173.4 keV ec end- 
point) followed by the decay of 85m Rb to 85 Rb (gamma 
of 514keV cc ) with 1.46 /is half-life and a branching ratio 
of 0.434%. Figure 15 shows the waveform of a 85 Kr can- 
didate event. The natural Kr concentration, derived with 
this method, is compatible with the one inferred from the 
data/Monte Carlo comparison. 
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FIG. 15: SI peaks of a candidate Kr event where the second 
light signal from the 7-ray is delayed by ~ 900 ns. 



The electronic recoil background in the signal re- 
gion consists of two contributions: leakage from the 
Gaussian shaped bulk of the electronic recoil back- 
ground distribution and a small contribution of anoma- 
lous (non-Gaussian) leakage. The first is determined 
by extrapolating the low-energy science data of the 
non-blinded log 10 (S2b/Sl) region. The prediction is 
(1.14±0. 48) events with the error being mainly due to 
the uncertainty in the definition of the electronic recoil 
discrimination level. 

Events with incomplete charge collection contribute to 
the anomalous background. "Gamma-X" events as de- 
fined in Sec. Ill E 5 have such topology, /^-decays from the 
dominating 85 Kr cannot create spatially-resolved multi- 
ple scatters due to the short electron mean free path, 
much less than the 3 mm position resolution. There- 
fore, 60 Co gamma calibration data was used to predict 
the anomalous background. The spatial distributions of 
anomalous events were found to be similar within 10% 
in the low energy region in Monte Carlo simulations of 
background and 60 Co calibration data. This was also 
confirmed by comparing spatial and energy distributions 
of potential leakage candidates, selected using the SI pat- 
tern likelihood parameter, in the 60 Co calibration data 
and low-energy non-blinded science data. Taking into 
account the different exposure of the 60 Co calibration 
compared to the science data, and the fraction of the 
background due to gamma interactions that can undergo 
multiple scatters, the anomalous leakage into the signal 
region was estimated to be (0.56±g;27) events. The er- 
ror includes the anomalous leakage discrepancy in the 
calibration and background distributions and the uncer- 
tainty in the estimation of the 85 Kr concentration in the 
LXe. 



3. Total background expectation and side-band test 

Taking into account all background sources mentioned 
above, the total background prediction for the bench- 
mark WIMP search region, 100.9 live days exposure and 
48 kg fiducial mass is (1.8 ± 0.6) events. Table I summa- 
rizes the background contributions. Alpha decays in the 
LXe do not contribute to the background as they have 
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TABLE I: Contributions to the total background prediction 
for the benchmark WIMP search region, 100.9 live days ex- 
posure and 48 kg fiducial mass. 

Source Expected background 

ER Gaussian leakage (1.14±0.48) events 
ER anomalous leakage (0.561q 27) events 
NR from neutrons (0.11±g g|) events 
Total (1.8±0.6) events 



Sec. HIE 5) as a function of nuclear recoil energy E nr . 
Based on the comparison of the observed event distribu- 
tion with the background-like distribution from Monte 
Carlo simulations, the background-only hypothesis was 
tested using a Profile Likelihood method [22]. The p- 
value for this hypothesis, which is the probability that 
the outcome of a hypothetical XENON100 experiment 
results in an observation equal or less background-like 
than the observed one, was found to be 31%. 



energies of a few MeV, far above the energy region of 
interest for the dark matter search. 

Before analyzing the science data in the region of in- 
terest, the side-band region (30—130) PE was unblinded. 
The background predictions for the 99.5%, 99.75% and 
99.9% discrimination levels were 14.3^gg, 8.2j^'g and 
5.61°'q events while 10, 8 and 7 events were observed, 
respectively. The observed and predicted numbers of 
events are compatible within the errors confirming the 
validity of the background predictions. 

Si [PE] 
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FIG. 16: The electronic/nuclear recoil discrimination param- 
eter as a function of nuclear recoil energy. Blue points in- 
dicate the observed event distribution after all cuts. Cyan 
points show the nuclear recoil distribution as measured with 
an 241 AmBe neutron source. The dashed-lines represent 
the energy window 8.4 — 44.6keV nr and the threshold cut 
S2 > 300 PE. The dotted-lines indicate the 99.75% rejection 
line from above and the 3a contour of the NR distribution 
from below. These are not used in the Profile Likelihood 
analysis, but are used to define the benchmark WIMP region. 
Three events with energies of 12.1, 30.2 and 34.6 keV llr fall 
into the benchmark region (red circles). This plot shows the 
same event distribution as Fig. 3 in [4]. 



G. Observed event population 

Figure 16 shows the final distribution of events in the 
electronic/nuclear recoil discrimination parameter (see 
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FIG. 17: Waveforms of the three events in the benchmark 
region that pass all cuts, including the ones defined post- 
unblinding. The SI peaks are labelled by open triangles, while 
the S2 peaks are labelled by full triangles. The insets show 
an expanded view around the SI region. Periodic electronic 
noise is visible in the top event, while single electron pulses 
after the major S2 at the trigger position (220 /is) are visible 
in the bottom event. 

In the pre-defined benchmark WIMP search region, 
three events remain after all cuts. These events with en- 
ergies of 12.1, 30.2 and 34.6 keV nr , as derived from their 
SI signal, are marked by red circles in Fig. 16. The events 
are homogeneously distributed inside the fiducial volume, 
as can be seen in Fig. 4 of reference [4] . Their waveforms 
are shown in Fig. 17. Though one of the events shows 
presence of periodic electronic noise in the waveform (top 
panel), the SI signal has a higher amplitude and is not 
affected by it. The other two events have no noise present 
in their waveform (middle and bottom panels). 

Within the background expectation of (1.8 ± 0.6) 
events, the 3 observed events do not constitute an ev- 
idence for dark matter as the probability of a Poisson 
process with an expectation value of 1.8 events to fluctu- 
ate up to 3 or more events is 28%. This is consistent with 
the p-value of the Profile Likelihood model. Therefore, 
exclusion limits on spin-independent elastic and inelastic 
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WIMP-nucleon scattering interaction were placed in [4] 
and [5]. 

IV. SUMMARY AND OUTLOOK 

A detailed description of the analysis of the 100.9 live 
days of science data acquired by XENON100 in 2010 is 
presented. The data set, the selection criteria, the eval- 
uation of the cut acceptances and the background esti- 
mates have been described. This data had been used 
to derive results in terms of spin-independent elastic 
[4] and inelastic [5] WIMP-nucleon interactions. How- 
ever, the methods presented here are general in nature 
and also apply to forthcoming XENON100 dark matter 
search data. Several improvements have been made on 
data quality cuts against electronic noise and position- 
dependent corrections of the measured signals using the 
significantly larger amount of both electronic and nuclear 



recoil calibration data. The most stringent limits to date 
on both spin-independent elastic and inelastic WIMP- 
nucleon cross sections demonstrated already the poten- 
tial of the XENON100 detector. A further refinement 
of the analysis, an improved detector performance and a 
greater exposure will increase the level of sensitivity of 
the XENON100 experiment. 
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